Draft version June 5, 2012 

Preprint typeset using 1^1^^ style emulatcapj v. 11/26/03 



(N 

O 
(N 

^: 

(N 



6 



>: 
o. 

On- 

o 



X: 



IMPLODING IGNITION WAVES. I. ONE-DIMENSIONAL ANALYSIS 

DoRON Kushnir\ Eli Livne^, and Eli Waxman^ 

Draft version June 5, 2012 

ABSTRACT 

We show that converging spherical and cylindrical shock waves may ignite a detonation wave in a 
combustible medium, provided the radius at which the shocks become strong exceeds a critical radius, 
i?crit- An approximate analytic expression for i?crit is derived for an ideal gas equation of state 
and a simple (power-law- Arrhenius) reaction law, and shown to reproduce the results of numerical 
solutions. For typical acetylene-air experiments we find -Rcrit ^ 100 /im (spherical) and i?crit ^ 1 mm 
(cylindrical). We suggest that the deflagration to detonation transition (DDT) observed in these 
systems may be due to converging shocks produced by the turbulent deflagration flow, which reaches 
sub (but near) sonic velocities on scales ^ i?crit ■ Our suggested mechanism differs from that proposed 
by Zel'dovich et al., in which a fine-tuned spatial gradient in the chemical induction time is required to 
be maintained within the turbulent deflagration flow. Our analysis may be readily extended to more 
complicated equations of state and reaction laws. An order of magnitude estimate of -Rcrit within a 
white dwarf at the pre-detonation conditions believed to lead to Type la supernova explosions is 0.1 
km, suggesting that our proposed mechanism may be relevant for DDT initiation in these systems. The 
relevance of our proposed ignition mechanism to DDT initiation may be tested by both experiments 
and numerical simulations. 

Subject headings: hydrodynamics - shock waves - supernovae: individual (la) 



1. INTRODUCTION 



The nature of the physical mechanisms driving a 
deflagration to detonation transition (DDT) is a ba- 
sic open question i n combust i on th e orv (see iLed 119771: 
Lee fc MoenI [l978l: IWihiamsl [19851 [Lewis fc Von Elbd 



1991 iHoflich et all [l995l : iHoflich fc Khokhlovl [1995 



Iwamoto et al.l[r999f ) are claimed to fit the observational 



1987t [Shepherd fc Led 11993 \K^dm^ for reviews) . 



DDT in the laboratory has been observed in sev- 
eral types of experiments in chambers contain- 
ing cold, unreacted, exothermic gases. These in- 
clude igniting the mixture with a spark (with or 
without obst acles along the walls of the chan- 
nels, _ see, e.g., lOppenheim et al.lll962t iLaderman et al.l 



1963; Urtiew & Oppenhei mI Il96a iP eraldi et al. 198L 
Teodorczyk 1995; Dorofce v et al.l]l996: Dorofcev 20oJ) 
or c reating a series of shock-flame interactions 
e.g.. iMarkstein et al.l [19641: IScarinci fc ThomasI [19901: 



Scar 



nci et al.[[1993[ : [Thomas et al.l[1997[ : [Thomas et"aLl 



20011 ). Experiments attempting to observe unconfined 
DDT (e.g., Knystautas et al. 1978) showed that a tran- 
sition to detonation induced by turbulent fiames in sys- 
tems without walls or obstacles is rather difficult. 

The most extreme case where a DDT may occur 
is in a Type la supernova (S NIa) e xplosio n, which is 
comp letely unconfined (see iHillebr andt fc Niemeved 
[2000 . for a review). One-dimensional detonation 
models of SNIa (I ArnettI 119691: [Hansen fc Wheeled [19691) 
do not reproduce the observed spectra. Both one- 
dimcnsional deflagration models ('Nomoto e t al.l [1976L 
[1984: Wo osley fc Weave r 1986) and delayc d-detonation 



Khokhlovl [19911: lYamaoka et al.[ 


1992[: iKhokhlov et al.l 


199.1 lArnett fc Livnd 1994allh 


; Wooslev fc Weaved 



^ Department of Particle Physics and Astrophysics, Weizmann 
Institute of Science, Rehovot 76100, Israel 

^ Racah Institute of Physics, Hebrew University, Jerusalem, Is- 
rael 

^ Similar to a model previously proposed by A. M. Khokhlov. 



data. 

DDTs observed in the laboratory are commonly 
characterized by turbulent flow flelds, initiated by 
the interaction of the laminar flame with walls, 
obstacles, or shocks, and by the initiation of 
a detonation wave by a sudden explosion ( "hot 
spot") . Two-dimensional ([Livnd[1993lJ: [Arnett fc Livnd 
[1994a; 'Khokhlov 1995"; 'Niemever fc Hillebrandt' '1995| 
Niemeyer et al. 1996;_Reinecke et al. 1999 ) and three- 
dimensional (iReinecke et al.l l2002ai bl: iGamezo et al.l 
[20031 [2004 [20050 SNIa flame simulations show that the 
gravity-induced Rayleigh-Taylor instability makes the 
flame turbulent, which allows extensive interpenetration 
of burnt and unburnt materials. The formation of hot 
spots in this case was not demonstrated computation- 
ally, presumably because of the small scales of the hot 
spots, which are not resolved in the simulations. 

The hot spots observed in the laboratory are believed 
to be pr oduced by a mechanism related to that proposed 
bv [Zel'dovich et al.l (|1970t ). in which a spontaneous re- 
action wave propagates through a reactive material with 
a spatial gradient in chemical induction time that leads 
to a deflagration prop agating at the Chapman- Jouguet 
(CJ) speed (see also Le e et al][1978[ ). However, it is not 
clear how the required flnely tuned spatial gradient is ob- 
tained. This "preconditioning" of the unburnt material 
appears to be artificial. Moreover, for both the labora- 
tory experiments and SNIa explosions, it is not clear how 
such a gradient can be maintained fixed in the turbulent 
deflagration flcld over the critic al length required for a 
successful ignition (see, however. lOran fc Gamezoll2007t ). 

Current studies of DDT use direct numerical simula- 
tions to investigate the possibility of the initiation of a 
detonation in turbulent deflagration flames. Due to the 
challenging computational requirements, calculations are 
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typically performed under simplifying approximations, 
and the initiation of detonation is typically not directly 
demonstrated. Rather, the turbulent flow field properties 
are studied and phenomenological arguments are used 
to argue f or the plausibility of the initiation of detona- 
tion (e.g.. IWooslev et alJ [20091 iPoludnenko et alJ 120111 : 
iWooslev et al.ll201lD . 

In this paper, we consider the ignition of a detona- 
tion wave by converging spherical and cylindrical shock 
waves, imploding onto an unburnt combustible medium. 
We limit our analysis to simple (ideal gas) equation of 
state and reaction laws, which allows us to investigate 
the various ignition modes and the flow behavior over 
a wide range of values of the relevant parameters. The 
model is described in detail in Section [2] The model 
equations are given in Section 12. 1[ dimensional analy- 
sis is used in Section to determine the characteristic 
length and timescales of the problem, and a summary of 
the dimensionless parameters is given in Section r2.3l The 
conditions under which a successful ignition is obtained 
are derived analytically in Section[21 and investigated nu- 
merically in Section [4] (the details of the numerical model 
are given in Appendix [B]). We find that the radius of the 
shock at fixed shock velocity should exceed a critical ra- 
dius in order to achieve a successful ignition. 

We note that ignition of detonation by converging 
shocks has been considered before, in studies of ig- 
nition by shock focusing thro ugh refl ection off con- 
cave walls (e.g . iBorisov et al.l [l 990: C han et all 119901: 
iBartenev et all 120001: iGelfand et al.i i2000l ) and by con- 
yerging toroidal shocks i n "pulse detonation engines" 
(iLi fc KailasanathI [200l iJackson et al.l 120031: iJacksonI 
I2005t ILI fc Kailasanatlj|2005[ ). However, these studies 
were mostly experimental, while analytical and numeri- 
cal investigations were limited to the study of the large- 
scale behavior of the flow in the experimental setup. In 
particular, a possible connection between the converg- 
ing flows and the generation of much smaller scale hot 
spots that may trigger detonation, as well as the hot spot 
properties required for detonation ignition, were not con- 
sidered. 

The analysis in this paper is one dimensional (limited 
to spherical or cylindrical symmetry). In a subsequent 
paper currently in preparation, we show that the evolu- 
tion of multidimensional perturbations during shock im- 
plosion does not suppress the ignition of a detonation 
wave. Although our analysis is limited to simple equa- 
tions of state and reaction laws, our predictions for the 
critical radius may be tested experimentally, since the 
properties of many explosive gases are well described by 
our simplified model. Moreover, our analysis may be 
readily extended to more complicated equations of state 
and reaction laws. Such extension for the case of pre- 
detonations in SNIa explosions, taking into consideration 
the effects of the large differences between the reaction 
lengths of the different burning stages in that case (see 
Section [S|), is presented in a third paper (in preparation). 

Our results are summarized and their possible impli- 
cations to DDTs are discussed in Section [S] 

2. THE MODEL 

In this section, we describe the model in detail. The 
model equations are given in Section 12.11 dimensional 
analysis is used in Section 12.21 to determine the char- 



acteristic length and timescales of the problem, and a 
summary of the dimensionless parameters is given in Sec- 
tion O 

2.1. Model Equations 

The reactive Euler equations, describing one- 
dimensional flow of an ideal gas that undergoes a single 
irreversible reaction, are 

{dt + udr) Inp -f r-^^-^^drir^-^u) = 0, 
{dt + udr)u + p^^drP = 0, 

{dt + udr)X = W, (1) 

where v — 1,2,3 for planar, cylindrical, and spherical 
symmetry, respectively. The total internal energy of the 
gas is given by 

P 



eip,V,X) 



XQ, 



(2) 



P(7 - 1) 

where p is the pressure, p is the density, u is the fluid 
velocity, A is the fraction of burnt material, 7 is the ratio 
of specific heats, and Q is the heat release per unit mass. 
We assume a reaction rate of the form 



(3) 



W = k[-^] (1 - A)' 



where p^, k, n, to, and pA are parameters characterizing 
the reaction. The argument of the exponential in Equa- 
tion ([2]) is the usual Arrhenius law for the temperature, 
— Ta/T, where Ta is the activation temperature, written 
in terms of density and pressure. 

In the planer case {v = 1), the reactive Euler equa- 
tions have solutions, the so called ZND waves (iZel'Dovichl 
ll940|:|Von Neumann 1942; Doring 1941), where a strong 
shock ignites the fuel and the burning proceeds to equi- 
librium in a reaction zone behind the shock, while the 
energy released continues to drive the shock. In the 
case where the flow at the end of the reaction zone 
is sonic, the detonation is called a CJ detonation, and 
the detonation propagates at the CJ velocity, given by 
Dc3 = \/2Q{'-f'^ — 1). We define the dimensionless pa- 
rameter 

PA PA 



PoDqj (7 + 1)PCJ 



(4) 



where pcj — PoD^j/ij + 1) is the CJ pressure (achieved 
at the end of the reaction zone), and po is the initial den- 
sity of the gas (without loss of generality, we may assume 
P-t, = Pq, where difference in their values is absorbed into 
K and into pa)- t measures how gradually the fuel burns 
in a detonation wave, as explained below. Let us define 
2^0.5(0.9) a-s the distance behind the shock front for which 
0.5 (0.9) of the fuel is burned (the reaction zone length 
is xi). Then, for small values of r the fuel burns slowly 
behind the shock (xo.g/xo.s — ^ 00), while as t 00 the 
profile of the burnt fraction behind the shock tends to 
a step function, in which a reactionless induction zone 
is terminated by a very thin "fire" in which all chemical 
energy is released (xo.g/xo.s — > 1). 

The initial conditions consist of a strong (spherical or 
cylindrical) shock wave, imploding from infinity, onto an 
unburnt stationary fluid with an initial density po- 
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2.2. Characteristic Length and Timcscales 

. . ,. Table 1. Some Values of o 

Let us nrst consider the pure nydrodynamic case (i.e.,0F 7 and u 

(5 = 0). In this case, the equations contain no character- 
istic length scales, and as the shock converges t he flow ap- 
proac hes Guderley's self-similar (SIS) solution ()Guderlevl 
Il942f ). This solution describes the flow for both the im- 
ploding stage (before the shock reaches the center) and 
the exploding stage (after the shock is reflected from the 
center as a finite Mach number shock). In what follows 
we analyze these two stages, as both are relevant for the 
ignition of a detonation wave. We denote the imploding 
(exploding) stage with a subscript i (e). 

In the SIS solution, the shock position, R, is given by 
the time relative to the time at which it reaches the cen- 
ter, io, as a power law 



S, io, Cs, Ae, fp,e, fp,e, AND fu,e AS FUNCTIONS 



R(x(t- tof 



(5) 



Equivalently, the shock velocity is related to the shock 
position through R oc R^ , where S = {a — l)/a. The 
value of a is determined by requiring the flow in the 
implosion stage to include a sonic point at some position 
fs{t) — CsR{t)- The resulting value of a is smaller than 
one {S is negative), such that the shock accelerates to 
an inflnite velocity as it converges to the center. Some 
values of a, 6, and as functions of 7 and v are given 
in Table [T] In what follows, we deflne i?cj as the radius 
at which R — —Da during the implosion, and define 

i = R/Rcj, 



R 



A,D. 



CJ 



/ R 



(6) 



where j = i or j = e, Ai — —1, and < < 1 (some 
values of are given in Table [1]) . At Rcj the internal 
energy generated by the hydrodynamic shock is compa- 
rable to Q (see Equation ([9])). 

The SIS solution and the Rankine-Hugoniot relations 
determine the post-shock density, pressure, and (the lab- 
oratory frame) velocity to be 

_f 7+1 
Ps.j-Jp,j ^ _ ^Pa, 

2 -2 

Ps,j = /pj , -, PoR 



' fu 



'^^7+1' 
2 

'7 + 1 



R, 



(7) 



fp,i — fu,', 



1, 



where R is the shock velocity, /p.i 
/p,e > 1, fp,e > 1, and < fu,e < 1 (some values of /p,e, 
/p,e, and fu,e are given in Table [1]). 

We turn now to analyzing the Q > case. Let us in- 
troduce a hydrodynamical time scale, t^j, and an energy 
per unit mass scale, Q/j.j, characterizing the pure hydro- 
dynamic (Q = 0) flow, and a timescale characterizing the 
chemical energy release, tqj. As long as 

_Q 

the chemical energy generation term in the third equa- 
tion of Equations ([T]) may be neglected, such that Equa- 
tions ([T|) reduce to the pure hydrodynamic Euler equa- 
tions, i.e., describing the pure hydrodynamic case. We 
assume that initially t^^i <C ig,i and that Equation ([8]) 



■inm{tqj,th,j) < Qh,: 



(8) 
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1.25 


4/3 


5/3 


1/ = 2 a 


0.85 


0.84 


0.82 





-0.17 


-0.19 


-0.23 
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fi 1 fi 

U. lU 


n 99 


1 .uu 


/- 


i.iO 


l.iD 


1 9n 


Ae 


0.21 


0.29 


0.58 


fp^e 


32.5 


18.0 


5.7 




822 


272 


29.5 


fu,e 


0.17 


0.18 


0.21 


u = 3 a 


0.74 


0.73 


0.69 


S 


-0.34 


-0.37 


-0.45 


€0 


0.31 


0.47 


1.00 


Cs 


1.11 


1.12 


1.15 


Ae 


0.21 


0.30 


0.64 


fp,e 


91 


38 


8.1 


fp,e 


1.65 X 10^ 


411 


30.6 


fu,e 


0.22 


0.24 


0.29 



holds through the entire flow, such that the flow ap- 
proaches the SIS solution. We define Qh,j as the post- 
shock internal energy excluding chemical energy. 



Qh,j 



We further define 

th. 



P 



2R-' 



X7-I) 
f, 



/p.. (7 + 1)2 



P,3 



A,D. 



CJ 



2S 



aR _ fh,ja Rcj 



A, Dr 



1-5 



(9) 



(10) 



where the dimensionless fh, 



^3 ^CJ 
is determined as follows. At 
the implosion stage, the available hydrodynamical time 
for a fiuid element to release chemical energy sufficient 
for affecting the fiow is the crossing time between the 
shock front and the sonic point (the flow behind the sonic 
point cannot affect the flow in front of it), which can be 

estimated as ~ -(Cs - l)i?/us,j = -(7 + l)(Cs - l)i?/2/^. 
This estimate leads to the deflnition 

/v = -(7+l)(Cs-l)/2a. (11) 

At the explosion stage the entire flow behind the shock 
is subsonic, and the hydrodynamical timescale may be 
estimated simply as the time passed from the shock's 
reflection, i.e., fh,e = 1- 
The post-shock chemical energy release timescale is 

" Up " 



t„ 



IdX 



Po 



1 



7+1 



jp.j J 



where 



1) 



(12) 



(13) 



ry = r(7-f 1)V2(7- 
If the inequality of Equation ([5]) is violated during the 
implosion, deviations from the pure hydrodynamical so- 
lution, e.g., ignition of a detonation wave, may arise. If 
Equation dH) holds throughout the implosion, no signif- 
icant deviation from the hydrodynamic solution arises, 
and a detonation wave is not ignited. In this case, the 
inequality of Equation ([S]) holds at the onset of explo- 
sion and may be violated as the shock explodes and Qt 
decreases. 
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2.3. Dimensionless Parameters 

The dimensional parameters defining the problem are 
Po, Rcj, Dcj, K, and pA- Based on dimensional argu- 
ments, these five-dimensional parameters may be used 
to construct three independent dimensional parameters 
(using which a dimensionless parameter may not be con- 
structed). Hence, the solutions are completely deter- 
mined (up to scaling of the measurement units of length, 
time, and mass) by a set of six dimensionless parame- 
ters: v, n, TO, and 7, which appear in the equations, 
and two additional parameters constructed from the five- 
dimensional parameters. We choose the latter two to be 
T (or 77, which is linear in t, see Equation (jl3p ) and 



61 = In 



where 



KaRcj 



n In 



4(7 - 1) 
7 + 1 



7 + 1 



-1/25 



2dln^o, (14) 



(15) 



Some values of ^0 & function of 7 and v are given in 
Table m 

The motivation for the choice of 9, which is essentially 
the ratio of i?cj to Dcj/k, is clarified in Section |3l We 
show in particular that ignition is obtained for values 
of 9 exceeding a minimum, or "critical", value 9c, and 
derive 9c as a function of the other five dimensionless 
parameters. The critical value of 9, 9c, corresponds, us- 
ing Equation to a critical value of i?cj- That is, a 
successful ignition requires Rcj (the radius at which the 
shock velocity reaches £>cj) to exceed a critical value. 
Using Equation (jG)), the lower limit for Rcj may be ex- 
pressed as a lower limit for the radius at which the shock 
reaches any chosen velocity. 

3. SIMPLE ANALYTIC ESTIMATES 

In this section, we present an approximate analytic 
analysis of the conditions under which a successful ig- 
nition is obtained. We define the ignition criterion, show 
that it implies that ignition is obtained for values of 9 
exceeding a minimum value 9c, and derive 9c as a func- 
tion of the dimensionless parameters. A detailed analysis 
of the flow behavior for different values of 9 is given in 
Appendices lA.ll and IA.2I for the implosion and explosion 
stages, respectively. In particular, the radii at which ig- 
nition is expected to be achieved are estimated. 

We define our successful ignition criterion as achieving 
chemical energy release that exceeds the hydrodynamical 
thermal energy generation at some shock radius. This 
requirement may be written as 



Q 



(16) 



where ft.j and fqj are order unity functions of the di- 
mensionless parameters (to be calibrated from numeri- 
cal simulations. Section |4]), and the hydrodynamic time, 
th,j, the combustion time, tq_j, and the hydrodynamic 
energy release, Qhj, were defined in Section [2?2] Note 
that thj_^ tqj for both ^ and ^ 00 (see Equa- 
tions (flOl) and (fT2l) ). and that the condition tq^ < th,j, 
which may be satisfied within some ^ range, is neither 
necessary nor sufficient for satisfying the ignition crite- 
rion of Equation (jl6p . 



During explosion, the shock may propagate into mate- 
rial which underwent significant combustion, even if igni- 
tion was not achieved. For the explosion stage we there- 
fore require the criterion of Equation (jl6p to be satisfied 
at a radius where the refiected shock propagates into un- 
burnt material, i.e., at a radius for which tq = X/{dX/dt) 
at the upstream of the shock is larger than t^^c (see Sec- 
tion |X2]). 

The requirement of Equation (|16|) is a necessary re- 
quirement, since otherwise the chemical energy release 
does not affect the fiow (Equation ^ holds), but not 
necessarily a sufficient one. In Section |4l we demon- 
strate that the inequality of Equation (fT6|) is actually 
a sufficient condition by investigating numerically a few 
representative cases: n = 0, to = 1, 7 = 5/3, r in the 
range [0.05, 2] for both cylindrical {ly = 2) and spherical 
{v — 3) geometries (numerical results for \n — 1,to = 
1,7 = 1.25, r = 0.7, 9 ~ 9c,v = 2, 3}, values relevant for 
laboratory experiments, are given in Section [5]). We also 
calibrate in Section |4] the values of ft,j and fqj for the 
examined parameters. 

As noted above, th^j <^ tqj for both ^ and ^ 00. 
For ftjthj < tq j E quation (1161) may be written, using 
Equations ©, ^Ui), (HI]), dH]), and dH]), as 



> 



- (1 - 35) In 5,(0- (17) 

Jp,j 

The relation between tq_j and t^.j is given by 
l^f?^) = .9,(0-A,-2<51n(X 

For tqj < ftjthj, which may be obtained during im- 
plosion/explosion. Equation (jl6|) may be written, using 
Equations ^ and (flS)) . as 



(18) 



C > Co,j = fQ,jA 



\2 Ip-j 

' f ■ 



-1/25 



Co. 



(19) 



Thus, Equation (fT6|) is satisfied if (and only if) Equa- 
tion (|17l) holds for some ^ > . 

An illustration of the function gj{S,) is shown in Fig- 
ure[TJ gji^) diverges for both ^ — (as ~ —(1 — 3S) In^) 
and ^ — 7* cx) (as ~ vfpjfpj^j'^C^^^), and its minimum 
is obtained at 



-2(577 



1^2 



-1/25 



(20) 



Defining ^c.j = iiiax{^„ij, ^oj}, the inequality of Equa- 
tion (fTO)) may be written as a lower limit for 9 (see Fig- 
ure [1]), 



Defining 



^c.j is given by 



_ 1 -3(5 25 .-1 



-2(5 



^ i Co J, V > Vo,j 



(21) 



(22) 



(23) 




r '• 



I I 

Fig. 1. — Qualitative illustration of the function gj{^), and of 
the significance of $m,j; So,j; ^^<^ ^c,j (see Section[3ll. 



Successful ignition 




T 



Fig. 2. — Qualitative illustration of the critical curves d^ jij), 
Equation I I21II . t measures how gradually the fuel burns (see Equa- 
tion (|4]|), and d essentially measures the ratio of Rcji the radius 
at which the shock velocity reaches Dqj ( see Equation I©), to 
Dq^^Ik (see Equation II14I I). »?o,i) Equation II22I I. is related to tqj 

by r,oj = TO,, (7 + 1)2/2(7-1)' 



Unstable 
detonation 



The inequality of Equation ([21]) is equivalent to that 
of Equation ([T6|) . For the implosion case, this implies 
that ignition is obtained provided 9 exceeds the critical 
value Qc,i- For the explosion case, a more stringent re- 
quirement for Q may be obtained since, as noted above. 
Equation (fTB)) is required to be satisfied at ^ for which 
tq at the upstream of the shock is larger than t/i^g. As 
explained in detail in Section IA.2| although this require- 
ment may affect significantly the radius at which ignition 
is achieved, it does not modify the critical value of d re- 
quired for ignition, 0c,e- In what follows, we therefore 
adopt Equation (|2T|) as the successful ignition criterion. 

For dc.i < Oc,e {9c,i > (^c,e), the critical ignition is 
obtained in the implosion (explosion) stage {ignition by 
implosion (explosion)) and 9c — 9c.i {9c = dc,e)- Note 
that for high enough the ignition is always obtained in 
the implosion stage. Thus, if for some set of parameters 
the critical ignition is by explosion, then by increasing 
9 a transition to ignition by implosion is obtained for 
9 = 9c,i- An qualitative illustration of the critical curves 
is given in Figure [2] 

For large values of r (77) , we have ^0 > which implies 
that for 9 > 9c the flow is significantly affected by chem- 
ical energy release at > = Co, where 9 = g{£,d) — A. 
We expect ignition to take place near ^ = Equal- 
ity, £.d ^ — Co, is obtained at critical conditions, i.e., 
9 = 9c- Thus, for large values of r and 9 = 9c, ignition 
takes place near R — ^qRcj with th ~ tq and Qh ^ Q- 
For small values of t, critical ignition takes place near 
C = Cm Of r^^/^l*!, which implies that for small r igni- 
tion is expected at i? i?cj with Qh/Q r <C 1. For 
more details see Section \K\ 

The analysis becomes more complicated for high val- 
ues of r, as the detonation wave becomes unstable 
to o ne-dimension al and multidimensional perturbations 
('see lSharpelll997l and references therein). It is not clear 
that this instability inhibits successful ignitions, since it 
can lead to some other stable configurat ion (i.e., multidi- 
mens ional detonation cell structure; see lFickett fc Dav"isl 
I2OOII for a review). In order to bypass the question of 
stability in this one-dimensional analysis, we assume that 
if a detonation wave formed and has propagated over a 
sufficient distance then the instability will lead to a stable 



configuration. Since the cell size of cellular detonations 
is ten to one hundred times the reaction zone length (for 
both laboratory experiments and SNIa explosions; see, 
e.g., TLcc 1984; Ganiczo ct al. 199^), we assume that a 
"sufficient distance" is one hundred times of the reac- 
tion zone length. This assumption should be examined 
in a multidimensional analysis and is beyond the scope of 
this work. However, our results do not depend strongly 
on the value adopted for the "sufficient distance" . 

Let us summarize briefly the estimates presented in 
this section. For a choice of v, n, m, and 7, there exists a 
critical curve in the (r, 0)-plane, given by Equation (j2ip . 
For 9 > 9c we expect a successful ignition, while for 
9 < 9c we expect a failed ignition. If the critical igni- 
tion is by explosion, then by increasing 9 a transition to 
ignition by implosion is obtained for 9 = 9c,i- For high 
values of r, the detonation wave becomes unstable, but 
we still expect that a detonation wave may propagate 
for a sufficient distance, provided 9 > 9c- In the example 
illustrated in Figure [2] the critical curve includes both 
ignition by implosion and ignition be explosion parts. In 
general, depending on the values of the dimensionless 
parameters defining the problem and of fqj, ftj, the 
critical curve may completely determined by a single ig- 
nition mode. The critical value 9c implies, equivalently, a 
critical radius for i?cj (see Equation ^^)- Hence, for a 
successful ignition, there exists a minimal value for i?cj , 
which in turn implies a minimal value for the radius of 
the imploding shock wave at a given velocity. 

4. NUMERICAL SIMULATIONS 

In this section, we use numerical simulations to study 
the behavior of the flow at different parts of the (r, 9)- 
plane (the details of the numerical model are given in 
Appendix [BJ. We compare the numerical results with 
the simple estimates of Section |3l which also allows us 
to calibrate the values of ftj and fqj- The numeri- 
cal setup includes a spherical (cylindrical) piston that 
moves inward into a cold ideal gas. The flow equations 
are solved numerically, with resolution sufficient for fully 
resolving the reaction zone. Each simulation is defined 
by the set of six dimensionless parameters v, n, m, 7, 
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r, and 9. We focus on the = 3, n = 0, m = 1, and 
7 = 5/3 case, and investigate the behavior of the flow in 
four different regions of the (r, 0)-plane. We first demon- 
strate a successful ignition by impfosion in Section 14.11 
Wc then demonstrate in Section 14.21 other behaviors of 
the flow in the (r, 6')-plane: failed ignition due to low 9, 
SL successful ignition by explosion, and a failed ignition 
due to instability of the detonation wave (high r). We 
explain the analysis of the (r, 0)-plane and summarize 
our results in Section 14.31 We show that the analytic 
estimate of the 9c^i critical ignition curve, Equation (|2ip . 
provides a good approximation to the numerical curve 
(better than 0.4 units in 9 for 0.05 < r < 2) for ft,^ ~ 1.7 
and fg^i ~ 0.88. Similarly, the analytic estimate of 9c,e 
provides a good approximation to the numerical curve 
(up to 0.6 units in 9 for 0.05 < r < 0.8) for /f,e ^ 0.41 
and fq^e — 1-7. The four cases discussed in Sections 14.11 - 
14.21 are marked in Figure [H We also show the results for 
the 1/ = 2 case (keeping all other dimensionless parame- 
ters the same) in Figure [71 

In the remainder of this section, we define t = as the 
time at which the shock reaches the center in Guderley's 
solution. We further normalize the time to tcj, which 
is minus the time at which the imploding shock reaches 

4.1. A Successful Ignition by Implosion: t = 1 and 
9 = 6.3 

In this section we study a typical case of a successful 
ignition by implosion, obtained for r = 1 and 9 = 6.3 
(we calibrate 9c = 9ci = 5.85 ± 0.05 for t — I, see Sec- 
tion US]). For ft^ ~ 6.41 and fg^ ~ 0.88 Equation ^ 
is satisfied at the implosion stage at ^ < — 1-15 (in 
this case — Cd,i2, see Scction lAiTj) . We therefore expect 
deviations from Guderley's SIS solution as the shock con- 
verges below A normalized pressure (p/po^cj) ni^-P 
(as a function of t/tcj and r/i?cj), which shows the ig- 
nition and the detonation waves, is shown in Figure [H 
The solid line in this figure represents A = 0.5. The spa- 
tial profiles of the normalized density {p/ po), normalized 
pressure {p/ pqDqj), and A at several times (marked with 
dashed lines in Figure ^ are shown in Figure H) 

At time t/tcj = — 3 the converging shock position is 
^ ~ 2.15 > £,d, and the dynamics are still not affected 
by the small chemical energy release (see A curves in 
Figure S]). As a result, the profiles are consistent with 
Guderley's SIS profiles (shown in red). Note that since 
Guderley's SIS profiles include a sonic point (red circle), 
the region between this point and the piston is super- 
sonic, such that deviations of the numerical solutio n from 
the SIS solution are to be expected (Waxman & S hvartsI 
fl993l ) . Such deviations are seen for the density profiles. 

After the shock implodes further to ^ ~ 1.35 > 
{t/tc! = ^1-5), the chemical energy release begins to 
be important, as the burnt fraction increases behind 
the shock. The numerical profiles begin to show large 
deviations from Guderley's SIS profiles. Shortly after 
(t/tcj = —0.3) two detonation waves are formed: one is 
moving outward (at r/i?cj — 1-45) into nearly unburnt 
material (see the small values of A in the upstream of this 
wave) and the second is moving inward (at r/ Rcj — 0.2). 
The latter hits the center shortly before t = (since the 
detonation wave is faster than the SIS imploding wave) 
and emerges as a regular shock wave (since the fuel in 




Fig. 3. — Normalized pressure map (logj^g(p/poD^j)) as a func- 
tion of t/tcj and r/Rcj, for u = 3, n = 0, m = I, -y = 5/3, t = 1, 
and 9 = 6.3. The solid line represents A = 0.5. The four times 
shown in Figure [4] are marked with dashed lines. 
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Fig. 4. — Normalized density (p/po, blue), normalized pressure 
{p/ PqDqj, black), and A (brown) as functions of r/i?cj at several 
times (t/tcj = —3, —1.5, —0.3, 3), for i/ = 3, n = 0, m = 1, ■y = 
5/3, T = 1, and d = 6.3. Guderley's SIS profiles are also shown at 
t/^CJ = — 3 (red), along with their singular point (red circle). 



the upstream for this shock is already burnt). This regu- 
lar exploding shock is seen at r/Rcj — 2.1 at t/tcj = 3. 
Note that at this time the pressure tends to a constant 
value toward the center, while the density tends to zero. 
This means that the temperature diverges in the center, 
very similarly to Guderley's SIS solution for the reflected 
shock. At this time the outgoing detonation wave is al- 
ready at r/Rc] ~ 3.1. This detonation wave meets our 
criterion for a successful ignition (see Section H75)) . 

4.2. Qualitatively Different Behavior of the Flow in 
Different Parts of the (r, 9) -plane 

In this section we investigate the flow in three cases, 
which are different from the successful ignition by implo- 
sion presented in Section HTT] The first case is a typical 
case of a failed ignition due to low 9: t = 1 and 9 = 4.3 
(we calibrate 9c — 9c.i — 5.85 ± 0.05 for r = 1, see Sec- 
tion 14. 3p . The normalized pressure map for this case is 
shown in Figure [5] (left panel) . The material begins to 
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Fig. 5. — Normalized pressure maps (log^oCp/pofcj)) ^ func- 
tions of tjtQj and t/Rqj, for = 3, n = 0, m = 1, and 7 = 5/3. 
The color map is the same as in Figure |3] and the solid line rep- 
resents A = 0.5. Left panel: r = 1 and = 4.3. Middle panel: 
T = 0.2 and 6 = 0.3. Right panel: r = 2 and 9 = 9.3. 



burn when the shock is at r ~ i?cj> but aUhough the 
burnt fraction increases behind the shock, the chemical 
energy release cannot affect the hydrodynamics, as it is 
too small. As a result, no detonation wave is produced, 
and the converging shock propagates to the center with 
roughly the same speed as in Guderley's solution (as can 
be seen in Figure EJ the shock hits the center at f ~ 0). 
After the converging shock hits the center, it emerges as 
a regular shock wave and later on begins to decelerate. 
As a result, the pressure decreases, the reaction is sup- 
pressed and the ignition fails. In this case our criterion 
for a successful ignition is not met (see Section 14.31) . 

The second case is a typical case of a successful igni- 
tion by explosion with r < tq: t = 0.2 and 6 — 0.3 (we 
calibrate 6^ = 0c,e = -0.85 ± 0.05 and ec,^ = 0.65 ± 0.05 
for T = 0.2, see Section |473|) . A normalized pressure map 
for this case is shown in Figure [5] (middle panel) . We 
expect ignition to occur during the explosion stage, at 
a radius ^ ~ 7 beyond which significant combustion did 
not take place (for fg^ ~ 1.7 and ft^ ~ 1.7 we obtain 
fo.e — 3 .0, 6q ~ —0.71 < 9 and Ccq,e2 — 7.4, see Sec- 
tion [X]2|. Indeed, we find numerically that the chemical 
energy release during the implosion stage is too small 
for ignition (although the material is burning, see the 
solid line which represents A = 0.5), while the reflected 
shock releases sufficient chemical energy to ignite a suc- 
cessful detonation wave as it moves into unburnt material 
around ^ = 7. The resulting detonation wave meets our 
criterion for a successful ignition (see Section . 

The third case is a typical case of a failed ignition due 
to instability of the detonation wave (high r): t — 2 
and e = 9.3 (we calibrate Oc = 9.95 ± 0.05 for r = 2, 
see Section 14. 3p . A normalized pressure map for this 
case is shown in Figure [5] (right panel) . In this case al- 
though a detonation wave is formed, it fails after some 
time (t/tcj — 2.5). As explained in Section[31 this failure 
is attributed to inherent instability of detonation waves 
with high values of r. Before the failure of the detonation 
wave, it manages to propagate only ~ 45 reaction zone 
lengths, such that it does not meet our criterion for a suc- 



cessful ignition (we require propagation over 100 reaction 
zone lengths, see Section 143]) . Higher values of 6 lead to 
detonation waves that propagate longer distances, such 
that they meet our criterion for a successful ignition (see 
Section [43]). 

4.3. Numerical Analysis of the {t, 9) -plane 

In this section, we summarize the analysis of the (r, 6)- 
plane. Once a criterion for a successful ignition is chosen 
(see below), the methodology of analysis is as follows. 
For our numerical analysis, we represent the (t, 0)-plane 
as a discrete rectangular grid of (t, 9) points, separated 
by 0.1 in both the t- and the 6'-directions. For values 
of T in the range of 0.05 < r < 2.5 we searched for the 
minimal value of 0, Og, for which a successful ignition is 
achieved and for the maximal value of 9, 9f, for which 
ignition is not achieved. The critical value of 9 is then 
defined as {9s + 9f)/2. Since our grid resolution is 0.1, 
9c is determined to an accuracy of ±0.05, which implies 
a ~ 5% accuracy in i?cj (see Equation P^ ). If the 
critical ignition is by explosion (see below for the implo- 
sion/explosion criterion) we search for the minimal value 
of 9 for which the ignition is by implosion {9i) and the 
maximal value of 9 for which the ignition is by explosion 
{9e)- The critical value of 9 for ignition by implosion is 
then given by {9i + 9e)/2. 

The derived values of 9c,i and 9c,e are shown in Fig- 
ure [6] (x's mark the four cases discussed in Sections 14.11 
and l4.2|) . The analytic critical curves of 9c,i and 9c,e, ob- 
tained by using the analysis presented in Section [3l are 
shown in Figure [5] for ft.i = 1.7, /q,^ = 0.88, ft.e = 1-7, 
and fg e = 1.7. These curves provide a good approxima- 
tion to the numerically inferred 9c,i and 9c. e (ignition by 
implosion: up to 0.4 units in 9 for 0.05 < r < 2, ignition 
by explosion: up to 0.6 units in 9 for 0.05 < r < 0.8). 

The numeric and the analytic critical curves, 0c, i and 
9c.e, for v = 2 (keeping all other dimensionless param- 
eters the same) are shown in Figure [T] The analytic 
curves obtained for ft^i = 3, fq^i — 0.73, ft.e = 1-5, and 
fQ,e = 0.95 provide a good approximation to the nu- 
meric 9c^i and 9c,e curves (ignition by implosion: up to 
0.3 units in 9 for 0.05 < t < 2, ignition by explosion: up 
to 1.3 units in 9 for 0.05 < r < 2). 

The different functional dependence of 9cj on t above 
and below tqj , predicted by the analytic analysis (see 
Equations ((23]i and pT|) ) is well reproduced in the numer- 
ical simulations. The r dependence of ^c.ji the radius at 
which ignition is initiated for 9 = 9c. j (determined from 
the numerical simulations in a manner similar to that 
of the determination of ^^.^2 and ^oq,e2 in Sections 14.11 
and l4.2l respectively) is also consistent with the analytic 
predictions: S,c,j is roughly given by ^ p.j f or r > tqj, and 
increases consistently with Equation (jzOp for r < tqj. 

We define a numerical solution in which a successful 
ignition was achieved as a solution in which a detonation 
wave has propagated Mth = 100 reaction zone lengths 
(see Section[3]for the motivation for choosing A/th = 100; 
our results are not sensitive to the exact value of Afth, 
see below). The number of reaction zones propagated 
was defined as the integral of the distance over which 
the burning wave propagated as a detonation wave, di- 
vided by the instantaneous reaction zone width (defined 
as the distance between the positions where A = 0.1 and 
A = 0.9 in the simulation). In order to determine whether 
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Fig . 6 . — Numeric (dashed lines) and analytic (full lines, Equa- 
tion critical ignition curves in the (r, 6)-plane for spherical 
(i' = 3) shock, with a reaction law characterized by {n = 0, m = 1} 
and an ideal gas adiabatic index of 7 = 5/3. Blue and red lines 
show implosion and explosion critical curves, respectively. Error 
bars reflect the accuracy of our numerical determination of 6c, 
±0.05. The analytic curves are obtained from the analysis pre- 
sented in Section [S] using /t i = 1.7, fQ^i = 0.88, /t^e = 1.7, and 
/(5,e = 1.7. The circles denote roj. The four cases discussed in 
Sections 14. 11 - are marked with crosses. 
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Fig. 7. — Same as Figure[6]for cylindrical i^v = 2) shock. The an- 
alytic curves are obtained from the analysis presented in Section [S] 
using /t,i = 3, /Q,i = 0.73, ft,e = 1.5, and /g^e = 0.95. 



the burning wave at a given instant is an "acceptable" 
detonation wave, we used the foUowing criterion. We 
examined the hydrodynamical profiles at times when an 
outward propagating burning front exists (before it hits 
the inward moving piston), and measured the distance 
between the position where A = 0.5 and the outermost 
shock wave, WNum- We then inspected the upstream 
density and pressure of the outermost shock wave, and 
derived the reaction front structure (ZND profile) of a 
planer CJ detonation wave propagating into a medium 
with such density and pressure. From the ZND profile we 
obtained the distance behind the shock front for which 
A = 0.5, wzND- Since the detonation waves obtained 
in our simulations are not exactly at a steady state and 
are spherical (cylindrical), we expect some difference be- 
tween WNum and wzND- We chose to define the burning 
wave as a detonation wave if WNum/wzND < Wth = 3. 



For most of the cases the decision between failure and 
ignition was not sensitive to the values of Wth and Mth, 
since in the cases of success we obtained WNum/wzND — 1 
for many reaction lengths, and in the cases of failure we 
quickly obtained WNum/wzND ^ 1- However, for high 
values of r the derived value of 9c is somewhat sensitive 
to the values of Wth and Mth. 

In order to decide whether a successful ignition is 
achieved in a given simulation by implosion or by explo- 
sion, we examined the flow at the time for which the det- 
onation wave has propagated Mth reaction zone lengths. 
If a reflected (regular) shock wave exists and it hasn't 
reached the detonation wave up to this time, we define 
this case as an ignition by implosion, otherwise as igni- 
tion by explosion (typically, the discrimination can be 
made by inspecting the flow at much earlier times). 

5. SUMMARY AND DISCUSSION 

We have analyzed the conditions under which converg- 
ing spherical and cylindrical shock waves ignite a detona- 
tion wave in a combustible medium. We have analyzed 
both analytically and numerically strong shocks implod- 
ing onto a combustible medium with an ideal gas equa- 
tion of state and a simple (power-law- Arrhenius, Equa- 
tion ([3])) reaction law. Given the geometry of the prob- 
lem (v = 2,3 for cylindrical, spherical symmetry), the 
adiabatic index 7 of the gas, and the power-law indices 
({n,m}) determining the reaction rate, the problem is 
fully defined by two dimensionless parameters: r and 
(see Equations (|4]) and (fT4|) ). r measures how gradually 
the fuel burns, and essentially measures the ratio of 
i?cjj the radius at which the shock velocity reaches Dqj 
(see Equation ([6])), to Dcj/k (see Equation (fT4)) V We 
have analyzed the different ignition modes (Sections [3] 
and 1^ , and have shown that ignition is obtained for val- 
ues of 6 exceeding a critical value, 9c{t) (Equation ((2T|) . 
Figure [2|) , and derived approximate analytic expressions 
for the Oc{t) curves (Equation ([21]) ). The analytic ex- 
pressions provide a good approximation to the exact, 
numerically derived curves (see Figures [6] and [T]) . 

The minimum value of 9 required for ignition implies 
a minimum value of i?cj required for ignition. Equiva- 
lently, a minimum value of i?cj corresponds to a mini- 
mum value of the radius of the shock at any other chosen 
velocity (see Equation For example, defining the ra- 
dius at which the shock becomes strong as R2 , the radius 
at which the shock velocity is twice the speed of sound 
Co in the unburned material, we have 

«,.«c,,(|^)"" (24) 

(this esti mate cannot be used for much weaker shock 
waves; see iPonchaut et~aIll2006D . 

We have shown that ignition may be achieved either 
during implosion or explosion (following reflection at the 
center) of the shock wave. The ignition mode at the 
critical curve depends on the values of {n, m, z^, 7}, but 
for high enough values of 9 the ignition is always obtained 
during implosion. For large values of r the detonation 
wave becomes unstable but nevertheless can propagate 
for a sufficient distance to guarantee a successful ignition 
provided 9 > 9^- 

Although our analysis is limited to simple equations of 
state and reaction laws, our predictions for the critical 
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radius may be tested experimentally, since the proper- 
ties of many explosive gases are well described by our 
simplified model. Moreover, our analysis may be read- 
ily extended to more complicated equations of state and 
reaction laws (e.g., the case of pre-detonations in SNIa 
explosions). 

This mechanism for the initiation of a detonation wave 
may be responsible for DDTs. Since near sonic ve- 
locities are obtained in the turbulent flow fields typ- 
ical of DDT on scales much larger than R2 (see be- 
low), it is not unreasonable to suggest that the turbu- 
lent flow produces converging shocks leading to an igni- 
tion of a detonation wave. It is interesting to note that 
similar ideas regarding ignition by "hot spots" hav e al- 
ready bee n proposed by iBrinklev Jr fc LewisI ()1959 ') and 
lOppenheim et al.l ()1962[ ). but were abandoned without 
any quantitative analysis. In order to demonstrate the vi- 
ability of our proposed DDT mechanism, a detailed anal- 
ysis of the turbulent flow, which is beyond the scope of 
this work, is required. The existence of convergent flows 
around unburnt material needs to be demonstrated, as 
well as the insensitivity of the ignition to inhomogeneities 
of the flow fields (density, velocity, burned fraction, etc.). 

Let us estimate Rqj and R2 for a stoichiometric 
acetylene-air mix with density of 1.58 x 10~''gcm~'^ at 
roo m temperature, as typic al for laboratory experiments 
(see[0 ran fc rxamezol[2007l for review on relevant ex- 
periments, as well as for equation of state and reaction 
rate prameters). For this gas we have 7 = 1.25, r ~ 0.7, 
n = 1, K ~ 1.58 X 10**s-\ and Dcj ~ 1.87 x lO'^cms-i. 
Assuming ftj — fqj = 1 we find tq ~ 0.66 < r 
for the spherical case, implying = Co — 0-31 and 
£,c,e — Co,e — 0.21, from which we obtain 0c, i — 7.4 and 
9c.e — 0. We therefore find that critical ignition is by 
explosion {6c = 0c,e), as should be expected due to the 
low value of 7 — 1, which implies a large amplification 
of the density by the reflected shock, and to the finite 
value of n, which implies that the reaction rate is sen- 
sitive to the density (see Table [T] for the downstream 
density at the explosion stage). A numerical determi- 
nation of 9c, obtained by the method described in Sec- 
tion[4l yields 9c ~ 9c,e = 2.87 ± 0.05, which is in reason- 
able agreement with our simple estimates. Consequently, 
we find i?cj — 15/.tm and, using cq ~ 3.2 x lO'^cms"^, 
i?2 — 300 /im. Assuming /t ^ = fq j — 1 for the cylin- 



drical case, we have tq — 0.98 > r, ^c.i = im.i — 0.25, 
£.c,e = £.m,e ^ 0.35, 9c,i ^ 8.1, and 9c,e ^ l.l. Once 
again, critical ignition is by explosion (9c ~ 9c^e)- The 
numerically determined value, 9c = 9c. e — 1.77 ±0.05, is 
in a close agreement with our simple estimates, and we 
find i?cj — 5 fj-m and i?2 — 2 mm. 

The estimates given above are consistent with the rel- 
evant terrestrial experiments in two senses. First, the 
scale of i?2 is much smaller than the size of the experi- 
mental channels, typically tens of centimeters, over which 
near sonic velocities are obtain ed in the turbulent flow 
field typical of DDTs (see, e.g.. lOran fc GamezdIMl . 
Hence, it is not unreasonable to suggest that the result- 
ing turbulence may produce converging flows on a scale 
^ i?2. Second, i?c,i is too small to be fully resolved 
in three-dimensional simulations ( currently reach i ng ce ll 
sizes of only ~ 100 /.tm, see, e.g.. iGamezo et aLll2005() , 
which might explain why this ignition mechanism was 
not observed in simulations (recall that for t > tq and ig- 
nition by explosion we expect explosion at i? = ^o.e^cj)- 

Similar statements can be made regarding the pre- 
detonation phase of white-dwarfs in delayed-detonation 
scenarios of SNIa. For typical pre-detonation condi- 
tions, po = lO^'gcm"'^ and equal mass fraction of ^^C 
and ^^O nuclei, an order of magnitude estimate of i?2 
(corresponding to a shock velocity ~ 7 x lO^cms"^) is 
0.1 km, smaller than the typical scale of the nearly sonic 
turbulence. An order of magnitude estimate of i?cj is 
0.01 km, well below the resolution of numerical simula- 
tions (which current ly reach ~ 1 km resolution, see, e.g., 
Gamezo et al.l [2"005[ ). Note that in this case the deto- 
nation wave consists of three burning layers: a Carbon- 
burning stage, followed by an Oxygen-burning stage and 
a Silicon-burning stage. The above estimates relate to 
a successful ignition of a detonation wave with steady 
carbon and oxygen burning layers. A detailed analysis 
of the ignition of detonation under conditions relevant 
for SNIa explosions will be given in a subsequent paper, 
currently under preparation. 
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APPENDIX 
IGNITION RADII 

In Section [3l we presented an approximate analytic analysis of the conditions under which a successful ignition is 
obtained. In this section, a detailed analysis of the flow behavior for different values of 9 is given. In particular, the 
radii at which ignition is expected to be achieved are estimated. The implosion and explosion stages are analyzed in 
Sections lA.ll and IA.21 respectively. 

The Implosion Stage 

An illustration of the functions Fd.iiO = 9ii0 ~ and Fcq.iiO = 9ii0 ^ ~ 2(51n(C/Co.i) is shown in Figure [AEl 
Let us consider first the case 77 > 770, i. The horizontal lines representing different values of 9 may intersect each of the 
functions Fd^i{S,) and i^eq,i(C) up to two times. We denote the intersection points (in case they exist) as ^^ ^1 < ^^.^2 for 
Fd,i(C) and ^eq.ii < Ceq,i2 for i^eq,i(0- C > ^cqM Or ^ < ^cqM wc have ft^ith.i < tq,i and Equation is satisfied 
for 9 > Fd^i{£,). Since r] > ryo.i we have ^o,* > Cm,* and 9c,i = fd,i(Co,i) = ^'oq,i(Co,i)- For ^ > ^c,i the flrst intersection, 
i.e., the largest ^ intersection, is at C = Cd,i2, which implies that ignition is expected at ^ = Cd,i2- 

Similar considerations lead to the conclusion that ignition is achieved at ^ = £,d.i2 also for 77 < 770, i. Thus, for the 
implosion stage ignition is achieved with ft.ith.i < tq^i. Equality, ft.ith.i — tq,i, is obtained bX 9 = 9c and 77 > 770, i. As 
6* — > 00 we have ^^.42 — {9/7])^^/'^^ — > 00, and ignition is obtained for R 3> -Rcj- 
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The implosion stage, r/ > J7o,j 
Shock propagation direction 




The implosion stage, r] < jjo,j 
Shock propagation direction 



Cm,* 



Fig. A8. — Qualitative illustration of the functions = ~ (solid line) and foqjCC) = ~ ~ 2i51n(^/5oj) (dotted 

line), and of the significance of 5m, j, ?0,j, a^nd ^c,j for the implosion stage with rj > r)o ^ (left panel) and r] < r)o,i (right panel). Dashed 
lines represent different values of 8. At § = ?d we have 9 = Fj, which implies that the hydrodynamic energy release and the chemical 
energy release are similar, and at § = ^cq we have 9 = Fcq, which implies that the hydrodynamic and chemical timescales are similar (see 
Section |3j. 
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The explosion stage 

Successful ignition 



^q,e ft,e^h,e 
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Fig. A9. — Qualitative description of the flow behavior at different parts of the (t, S)-plane. We expect ignition near g = gig 



A qualitative description of the flow behavior at different parts of the (r, 0)-plane is given in the left panel of 
Figure EH Successful ignition is predicted at ^ = ,^ig — £^d.i2 for 9 > 9c,i- For large values of t (77), we obtain 
9c,i — V^Qi^ J ^iid critical ignition takes place near R — ^o,iRcj- For small values of t, critical ignition takes place 
near S^rn.i oc r^^/^l'^l » 1 and 9c,i oc Inr. In this case ignition is expected at i? ■fm.i^cj ~ t'^^^Dcj/k with 
Qh,,/Q r <Cl. 

The Explosion Stage 

As explained in Section [31 the exploding shock may propagate into burnt material, even if ignition is not achieved 
during implosion. Thus, for the explosion case we require Equation (|16p to be satisfied at a radius where the reflected 
shock propagates into unburnt material, i.e., into material for which tq in the upstream is larger than th,e- For 
simplicity, let us first assume that tq in the upstream of the shock is similar to tq in the downstream, (ig e). This is 
not a bad assumption since the reflected shock has a finite Mach number, so that the pre- and post-shock density 
and temperature are typically not very different. We later comment on modifications of our conclusions due to the 
inequality of the pre- and post-shock value of tq. 

Let us first consider the ?/ < 770. e case, which is illustrated on the right panel of Figure lAlOl and for which ^o,e < Cm.e 
and 9c.e — Fd,e{S,m,e)- For this case it is useful to define 69 = ^£;,e(Co,e) — -Peq,e(63,e)- For 9 > 9^ we expect that 
ignition would not be achieved at ^ = ^o,e, although Equation ([T51) is satisfied there, since the reflected shock is likely 
to propagate into burnt material at ^ = ^o,e- Instead, we expect ignition only at ^ = Coq,e2- For 9q > 9 > 9c,e, on the 
other hand, we expect ignition at ^d.ei- 

The r] > 7yo,e case, for which fo.e > ^m.i and 9c,e — ^d,e(Co,e) — ^cq.e(Co,e), is illustrated in the left panel of 
Figure lAlOl Arguments similar to those applied for the 9 > 9q and rj < 779, e case imply that for 9 > 9c,e ignition is 
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likely to be achieved near ^ = Coq,e2- For this case all ignitions are obtained with tq^e ft,eth,e- 

A comment is in place here regarding the assumed equality of tq in the upstream of the shock and tq e- Since the 
density and the temperature in the upstream of the shock are lower than the downstream values, tq in the upstream of 
the shock is larger than tq^e- This implies that the shock propagates onto unburnt material for some ^o,e < C < ^cq,e2, 
and therefore ignition is likely obtained in that range. As noted above, since the Mach number of the reflected shock 
is finite, the post and pre-shock conditions are typically not very different, so that ignition is likely to be achieved near 
Coq,e2- Furthermore, since ^o.e = foq.e2 for 0c,e, our estimates for the critical value of required for ignition are not 
affected by the inequality of tq in the upstream and tq^g. 

A qualitative description of the flow behavior at different parts of the (r, 6')-plane is given in the right panel of 
Figure EOl In this case ignitions are obtained with tq^e ^ ft,eth,e {fQ,eQh,e/Q > 1) for rj > ryo.e and for {9 > 9o,-q < 
ryo,e}- Otherwise, ignitions are obtained with ft,eth,e < tq ^- Above the critical curve for rj > ?7o,e, and above for 
V < Vo,ej we expect ignition near Ceq,e2- In the region between 6 = 9q and the critical curve for rj < 770, e we expect 
ignition at ^ = ^d,ei- For — > 00 we have ^cq,e2 — (/p,e^e^//p,e?/)~^^^'' — ^ 00, and ignition is obtained for R ;» i?cj- 
For large values of r (77), we obtain Oc^e — vfp,efpe^7^^ol^j and critical ignition takes place near R — ^o.e^cj- For 
small values of r, critical ignition takes place near ^jn,e cx r^^/^l"^' and, similarly to the implosion case, 6c,e oc Inr and 
ignition is expected at i? im,iRc3 ~ t^^^-Dcj/'* with Qhd/Q ~ r ^ 1. 

THE NUMERICAL MODEL 

We used the one-dimensional, Lagrangian version of the VULCAN code (for details, see iLivnil 1993al) . with a simple 
extension to include reaction rates of the form of Equation (jS]). The initial conditions used are zero velocity, constant 
density poi and small pressure po (see below) throughout the computation's region r < Lp. The initial mesh spacing 
was uniform with N cells, and we did not use any remapping or AMR technics. The outermost node acted as a piston 
with constant inward velocity, Vp, throughout the simulation. The initial pressure was chosen such that the initial 
sound speed was ~ O.lwp. This ensures that initially the ratio between the upstream pressure and the downstream 
pressure of the imploding shock wave is ~ 0.01. A further reduction of the initial pressure did not change significantly 
any of our results, but increased the computing time due to the Courant condition. We also suppressed to zero the 
reaction rate for p/ p < l.OOlpo/po in order to ensure that the material is completely unburnt until it is shock heated. 

The values of Lp, Vp, and po, which determine the dimensional scales of the problem, were chosen arbitrarily and 
were held fixed for all of our simulations. Given the values of u, n, m, and 7, we used the following consideration to 
set the values of k, pa, and Q in order to calculate the fiow for any choice of (t, 0). The velocity of the shock wave as 
it emerges from the piston is 

Hence 

fvVp^-Dcj(^^^ , (B2) 

and we have 

^CJ = -fvVpK-\ (B3) 

where Rqj = Lp/K (we discuss below how to choose the value of K). The CJ velocity sets the value of Q, and 
Equation (|4]) sets the values of pA- The value of k is given now by Equation p^ . 
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In order to determine if a successful ignition took place, we investigate the behavior of the outgoing waves. Our 
criterion for a successful ignition is that a detonation wave with WNum/wzND < wth (see Section |4|) has propagated 
over at least 100 reaction zone lengths. Since these waves ultimately hit the piston, we have to make sure that the 
outgoing waves either satisfy the criterion for a successful ignition or that WNum/ifZND diverges with time, before the 
piston is reach by the wave. Increasing the available distance for the outgoing wave (where at some point a decision 
regarding whether a successful ignition took place can be made) can be done by increasing the value of K. However, 
since we want to solve for the reaction zone behind a possible detonation wave with a given resolution, increasing 
K also increases N. The value of K was chosen to be the minimal value for which a decision regarding whether 
a successful ignition took place can be made, and the value of N was chosen such that the reaction zone was fully 
resolved (we always checked for convergence of our results, and in general we found that ~ 100 cells are sufficient for 
a converged calculation of the reaction zone). With the method presented here we were able to calculate all cases of 
interest with < 5 x 10^. 
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